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Abstract 

During  battery  discharge,  the  heat  generated  is  the  sum  of  the  Joule  (resistive)  and  enthalpic  (chemical)  heating  effects.  Conversely,  during 
battery  charging,  the  heat  generated  is  the  Joule  minus  the  enthalpic  heating.  If  the  conditions  are  carefully  selected,  one  can  observe  a  net  battery 
cooling  during  charging. 

However,  an  interesting  phenomenon  takes  place  during  overcharge.  Those  cells  designed  as  sealed  recombinant  systems  develop  significant 
heating.  Flooded  designs  do  not  exhibit  this  effect.  The  applied  electric  power  generates  energetic  reaction  products  as  a  consequence  of  the 
electrochemical  reactions.  This  is  an  energy  absorbing  process.  The  gasses  are  then  vented  into  the  environment.  Since  the  sealed  cells  undergo  a 
closed  recombination  cycle,  i.e.,  no  material  is  exchanged  with  the  environment,  the  rate  of  heat  generated  is  proportional  to  the  power  input  to  the 
cell.  Essentially,  the  cell  is  behaving  in  the  manner  of  a  resistor. 

In  this  connection,  the  thermal  runaway  phenomenon  that  has  been  often  observed  in  starved  electrolyte  cell  designs  raises  a  potential  problem 
in  battery  applications. 

It  is  not  efficient  to  design  around  the  worst  case  scenario,  i.e.,  anticipating  the  thermal  runaway  effect.  It  is  wiser  to  detect  its  onset  and  shut 
down  the  charging  process.  An  alternative  approach  is  to  develop  an  understanding  of  the  thermal  runaway  process  and,  perhaps,  develop  a  method 
for  eliminating  or  effectively  controlling  it. 

A  study  was  performed  in  an  attempt  to  model  the  thermal  runaway  effect.  In  short,  the  effect  appears  to  be  related  to  the  electrolyte  distribution 
in  the  separator.  This  suggests  that  modification  of  the  AGM  separator  properties  could  provide  a  means  for  better  controlling  the  thermal  runaway 
failure  mode. 
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1.  Introduction 

Thermal  runaway  is  a  prevalent  failure  mode  appearing  in 
recombinant  lead-acid  batteries  [1].  On  occasion,  VRLA  bat¬ 
teries  are  removed  from  service  because  the  battery  gets  quite 
hot  to  the  touch  during  constant  voltage  charging.  This  heating 
has  often  been  attributed  to  the  exothermic  recombination  of 
oxygen  at  negative  electrodes.  Under  normal  conditions,  con¬ 
stant  voltage  charging  of  lead-acid  batteries  shows  a  decrease 
in  current  approaching  an  asymptotic  limit.  In  the  case  of  the 
thermal  runaway,  the  current  can  rise  to  the  power  supply’s  cur¬ 
rent  limit.  Fig.  1  shows  a  normal  charging  curve  for  a  VRLA 
battery.  The  early  constant  current  portion  is  the  consequence 
of  the  power  limitation  of  the  power  supply  used  to  charge  this 
battery.  The  declining  portion  of  the  curve  is  the  consequence 
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of  the  increasing  battery  voltage  resulting  from  the  increasing 
state  of  charge  of  the  battery.  The  current  eventually  reaches  an 
asymptotic  limit  where  small  internal  parasitic  processes  sustain 
a  very  low  current  flow. 

In  contrast,  Fig.  2  shows  the  battery’s  response  during  thermal 
runaway  observed  while  the  battery  was  being  charged  under 
normal  constant  voltage  charging  conditions. 

When  this  failure  mode  appears,  the  effect  is  very  evident 
and  is  easily  identified.  An  examination  of  this  failure  mode  has 
value  in  order  to  understand  why  it  occurs  and  by  doing  that,  to 
determine  ways  of  minimizing  its  occurrence. 

2.  Experimental  note 

The  author  is  aware  of  two  methods  for  inducing  the  ther¬ 
mal  runaway  of  a  VRLA  battery  in  the  laboratory.  Application 
of  these  methods  provides  a  means  of  accelerating  the  thermal 
runaway  effect  for  the  purpose  of  study.  The  first  and  quickest 
way  to  induce  the  effect  is  to  charge  a  VRLA  battery  well  above 


978 


H.A.  Catherine)  /  Journal  of  Power  Sources  158  (2006)  977-986 


Normal  CV  Charging  Curve 
For  Hawker  Armasafe  Battery 


Seconds 

Fig.  1.  Normal  constant  voltage  charging  curve  for  battery. 


its  recommended  charging  voltage.  For  example,  a  VRLA  bat¬ 
tery  designed  for  a  nominal  12  V  operation  could  be  charged 
at  16  V  to  induce  the  effect.  The  second  method  takes  a  little 
longer  to  bring  about  the  effect.  Specifically,  a  VRLA  battery 
is  cycled  using  a  constant  current  during  both  charge  and  dis¬ 
charge.  The  charging  takes  place  at  about  a  C/2  rate  until  the 
temperature  at  the  negative  terminal  rises  to  a  maximum  of 
40  °C.  The  charging  rate  is  determined  by  noting  that  C  is  the 
capacity  in  Ah  and  the  denominator  is  in  hours.  The  subsequent 
discharge  can  be  carried  out  a  C/10  rate  to  an  appropriate  cut¬ 
off  voltage  (i.e.,  10.5  V  for  a  12  V  nominal  battery).  Repeating 
this  for  50  cycles  and  then  initiating  a  standard  constant  volt¬ 
age  charge  to  the  battery  from  a  discharged  state  has  resulted 
in  a  thermal  runaway.  In  the  study  here,  this  was  the  preferred 
method.  The  repetitive  cycling  serves  to  precondition  the  bat¬ 
tery  to  make  it  susceptible  to  thermal  runaway.  The  desired 
objective  was  to  demonstrate  the  effect  as  it  would  take  place 
under  “normal”  operating  conditions.  This  methodology  may 
most  closely  represent  what  appears  when  the  thermal  runaway 
failure  takes  place  under  normal  applications.  Of  course,  the 
usual  caveats  must  be  kept  in  mind  whenever  one  uses  acceler¬ 
ated  testing  methods  to  induce  a  failure  mode.  Also,  the  methods 
described  above  have  not  been  optimized  to  induce  any  particular 
result.  In  the  tests  performed,  the  rated  capacity  of  the  batteries 
was  nominally  120  Ah.  The  discussion  that  follows  attempts  to 
give  an  interpretation  to  the  variety  of  behaviors  observed  dur¬ 
ing  thermal  runaway  as  observed  using  the  induction  procedure 
described. 
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Fig.  2.  Current  and  temperature  response  of  a  VRLA  battery  in  thermal  runaway. 


3.  Discussion 

The  phenomenon  of  thermal  runaway  was  discussed  by 
Pavlov  as  a  consequence  of  the  closed  oxygen  cycle  in  valve 
regulated  lead-acid  batteries  [2] .  Essentially,  during  the  charg¬ 
ing  of  a  VRLA  battery  and  starting  at  about  70%  state  of  charge, 
oxygen  begins  to  evolve  at  the  positive  electrode  at  very  low 
rates.  The  design  of  the  battery  provides  an  efficient  means  for 
removing  the  oxygen  by  recombination  at  the  negative  electrode. 
This  charging  process  normally  takes  place  at  constant  applied 
voltage.  The  observed  charging  current  decreases  exponentially 
and  approaches  a  low  steady  state  value.  This  low  limiting  cur¬ 
rent  is  indicative  of  the  completion  of  the  charging.  The  oxygen 
recombination  process  continues  at  a  very  low  rate  at  the  applied 
charging  voltage.  This  is  the  normal  behavior  observed  during 
the  battery  charging. 

Pavlov  presented  a  critical  insight  into  the  thermal  runaway 
effect.  The  interpretation  is  based  on  an  observation  made  by 
Pavlov  and  Monahov  [3].  Essentially,  the  normally  operating 
VRLA  battery  is  a  closed  system.  In  the  closed  system,  the  rate 
of  the  oxygen  recombination  at  the  negative  electrode  is  deter¬ 
mined  by  the  rate  of  oxygen  evolution  at  the  positive  electrode. 
The  oxygen  generation  rate  is  dependent  on  the  potential  and 
the  temperature  at  the  interface  of  the  positive  electrode.  Repro¬ 
duced  here  is  a  Tafel  plot  showing  the  critical  dependencies  of 
the  electrochemical  oxygen  evolution  on  the  PbC>2  layer  on  the 
positive  electrode  at  various  measured  cell  temperatures.  The 
experimental  data  is  shown  in  Fig.  3. 

The  influence  of  temperature  can  be  seen  by  holding  the  elec¬ 
trode  potential  constant  and  observing  the  rate  of  the  oxygen 
evolution  reaction  which  is  proportional  to  the  log  of  the  mea¬ 
sured  current  density.  As  shown  in  the  figure,  by  holding  the 
electrode  potential  equal  to  1.3  V  versus  a  Hg/HgSCL  reference 
electrode,  the  currents  at  the  different  temperatures  can  be  iden¬ 
tified.  For  example,  at  42  °C  the  current  density  passing  through 
the  electrode  is  0.26  A  cm-2  as  compared  to  1.072  A  cm-2  at 
52  °C.  This  is  showing  that  a  10  °C  temperature  increase  is  caus¬ 
ing  a  four-fold  increase  in  the  rate  of  oxygen  evolution.  This  is 
a  characteristic  Arrhenius  temperature  dependence  on  rate  pro¬ 
cesses: 

rat  e  =  ke~AH/RT  (1) 


Rate  of  Oxygen  Evolution 


Fig.  3.  Tafel  dependencies  of  the  oxygen  electrode  at  different  temperatures. 
The  potential  of  the  oxygen  electrode  is  measured  vs.  a  Hg/Hg2S04  reference 
electrode. 
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In  this  case,  AH  is  the  activation  energy  of  the  reaction,  R  the 
gas  constant  and  T  is  the  absolute  temperature.  The  approach 
to  the  mechanistic  interpretation  discussed  below  is  based  on 
observing  the  exponential  dependence  of  the  oxygen  evolution 
reaction  on  temperature  while  holding  the  electrode  potential 
constant. 

3.1.  Heat  management  discussion 

To  begin  to  analyze  this  complex  process,  it  is  helpful  to 
consider  the  situation  where  a  fixed  current  acts  to  generate  oxy¬ 
gen  that  recombines  at  the  negative  electrode.  A  sealed  system 
is  assumed.  The  heat  thus  generated  establishes  a  steady  state 
with  the  heat  removal  process.  To  simplify  the  description  of  the 
modeled  process,  it  is  convenient  to  assume  that  the  charging 
component  associated  with  energy  storage  is  complete  so  that 
the  heat  production  is  represented  as  coming  from  Joule  heating 
and  the  enthalpy  of  reaction  of  the  closed  cycle  oxygen  recom¬ 
bination.  The  rate  of  heat  removal  is  given  in  its  usual  form 
as  being  proportional  to  the  temperature  gradient  between  the 
source  and  the  environment 

rate  of  Joule  heating  +  enthalpic  heating  =  rate  of  heat  removal 


Although  this  observation  is  nothing  more  than  a  requirement 
placed  on  the  system  by  the  energy  conservation  law,  a  few  addi¬ 
tional  comments  are  in  order  because  this  conclusion  stands  in 
opposition  to  some  widely  held  views  on  the  subject.  Specifi¬ 
cally,  it  is  a  commonly  held  view  that  the  oxygen  recombination 
step  is  the  heat  generating  mechanism  taking  place  during  over¬ 
charge  in  a  sealed  battery  system.  It  is  also  worth  noting  that 
various  authors  have  taken  guarded  positions  on  this  admitting 
only  that  a  correlation  exists  between  the  oxygen  recombination 
process  and  the  heat  generation  effect  [4-6].  The  “leap  of  faith” 
takes  place  as  the  result  of  the  well-established  observation  that 
the  reaction  of  oxygen  with  metals  is  known  to  be  exothermic. 

To  support  the  position  taken  in  this  exposition,  consider  the 
following  reactions  as  representing  a  closed  cycle  process.  (One 
can  add  more  complexity  if  desired  but  this  keeps  it  simple.) 

During  overcharge: 

•  At  the  negative  electrode 
2Pb  +  02  ->  2PbO, 

PbO  +  2H+  +  2e“  -*  H20  +  Pb  (2X) 

•  At  the  positive  electrode 


I2R  +  keT  =m(T  -T0)  (2) 

In  this  case,  R  is  the  resistance,  I  the  current,  T  is  temperature 
and  m  is  the  heat  transfer  coefficient.  It  is  to  be  noted  that  the 
magnitude  of  the  current  which  is  proportional  to  the  rate  of 
oxygen  production  is  in  turn  proportional  to  the  rate  of  oxygen 
recombination.  That  is 


/  =  K  e7 

(3) 

then,  substituting: 

(K  er)2  +keT  =  m(T  -  T0) 

(4) 

or 

k\  e2r  +  k2  er  =  m(T  —  7q) 

(5) 

The  first  term  on  the  left-hand  side  (LHS)  of  the  equation  can 
grow  more  rapidly  than  the  second  term.  However,  the  coeffi¬ 
cient  will  determine  the  magnitudes  of  each  quantity.  Further,  the 
value  of  R  may  well  be  increasing  in  magnitude  as  the  evolved 
gas  displaces  the  electrolyte  in  the  separator. 

At  this  point,  it  is  instructive  to  ask  which  term  on  the  left- 
hand  side  of  the  previous  equation  is  larger.  That  is,  is  it  the  Joule 
heating  or  the  chemical  enthalpy  that  predominates?  A  prevalent 
view  is  that  the  recombination  of  oxygen  is  said  to  be  the  major 
source  of  the  heat  generation  during  overcharge. 

It  is  important  to  emphasize  that  a  position  will  be  taken  here 
that  is  contrary  to  that  viewpoint.  That  is,  in  the  closed  oxygen 
cycle,  the  net  enthalpy  of  the  reaction  is  zero.  This  means  that  the 
closed  oxygen  cycle  contributes  no  net  heat  to  the  overcharge 
process.  It  is  to  be  understood  that  the  meaning  of  a  “closed 
cycle”  infers  that  no  material  leaves  the  battery  and  that  the 
reactants  and  products  return  to  their  original  state.  In  that  case, 
the  enthalpy  of  the  cyclic  process  is  necessarily  zero. 


2H20  -*  02  +  4H+  +  4e“ 


Then  this  sequence  of  processes  can  be  represented  as  a 
Born-Haber  cycle.  This  sequence  of  reactions  (Fig.  4)  can  be 
used  as  an  example. 

This  cycle  illustrates  a  sequence  of  the  states  of  matter  cor¬ 
responding  to  the  chemical  changes  occurring  during  the  cycle. 
To  illustrate  the  application  of  this  cycle,  the  following  exercise 
is  presented.  AH  is  the  enthalpy  of  each  particular  reaction 

A  Hi  +  A  H2  +  AH3  +  AH4  =  0 
(the  first  law  of  thermodynamics) 

AH\=0  (no  change  of  state  takes  place),  AH2  =  +  571.88  kJ 
mol-1  O2  [2],  AH3  =  —  (2)(195.2)  kJ  mol-1  Pb  (standard  heat  of 
formation  of  PbO),  AH4  =  quantity  to  be  determined. 

Then 

A H4  =  [0  -  571.7  +  (2)(195.2)]/2kJmol_1  PbO 
=  90.65  kJ  mol-1  PbO 


2Pb 

1  i 

AH! 


2Pb 


AH3 

02  - *  2PbO 


AH 


2 


-4H+ 
-  4e' 


+  4H+ 
+  4e' 


AH4 

+  2H20  * - 


Fig.  4.  Born-Haber  cycle  illustrating  the  essential  parts  of  the  closed  oxygen 
cycle. 
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for  the  reaction: 

PbO  +  2H+  +  2e~  ->  H20  +  Pb 

The  point  is  that  the  sum  of  all  the  enthalpies  associated  with 
each  of  the  participating  reactions  occurring  in  a  closed  cycle 
must  equal  zero.  Otherwise  one  can  create  or  destroy  energy 
depending  on  the  direction  one  moves  around  the  cycle. 

As  applied  to  the  closed  oxygen  cycle,  the  claim  that  no  net 
heat  is  generated  is  a  theoretically  based  argument.  No  empirical 
data  is  presented  here  to  support  it.  It  is  noteworthy  that  the 
position  taken  in  this  discussion  involves  a  different  paradigm. 
That  is,  either  energy  is  conserved  or  it  is  not  conserved.  The 
position  taken  is  based  on  first  principles  and  not  on  experimental 
data.  There  is  an  interesting  philosophical  issue  involved. 

As  a  summary,  the  exothermic  heat  of  recombination  of  oxy¬ 
gen  at  the  negative  electrode  is  acknowledged.  However,  all  of 
the  reactions  taking  place  in  the  electrode  pack  have  to  be  con¬ 
sidered.  The  sum  of  the  heats  of  reaction  of  each  of  the  steps 
leading  to  the  recombination  and  regeneration  the  reactants  is 
the  enthalpic  component.  In  that  case,  the  closed  cycle  neces¬ 
sarily  requires  that  no  net  heat  can  be  generated.  If  this  is  taken 
to  be  the  case,  then  the  only  source  of  heat  produced  during  the 
overcharge  is  the  Joule  heating  contribution. 

3.2.  Examination  of  non-linear  effects 

For  the  discussion  that  follows,  it  is  only  necessary  to  observe 
that  the  heat  production  rate  is  an  exponential  function  and  the 
heat  removal  process  is  a  linear  function. 

Taking  the  Joule  heating  as  the  only  source  of  energy  pro¬ 
ducing  heat,  the  Joule  component  can  be  given  as 


/Foe  T  /  ^polarization  —  m(T  To) 

(6) 

where  V0c  is  the  open  circuit  battery  voltage.  Also: 

/Foc  ^  /  ^polarization  dnd  I  =  K  e 

(7) 

and  substituting: 

(Kek"T)V0C  =m(T  -To) 

(8) 

or 

ki  ek"T  =  m(T  -  To) 

(9) 

It  is  the  behavior  of  the  mathematical  functions  in  the  steady 
state  relationship  that  is  instructive.  Fig.  4  shows  a  functional 
response  for  the  purpose  of  showing  the  steady  state  condition: 

k  ek" T  =  m(T  -  T0) 

00) 

Fig.  5  plots  the  exponential  rate  of  heat  generation  and  the 
linear  rate  of  heat  removal.  In  this  plot,  m  was  assigned  a  value 
of  5,  A:  and  kf  the  value  of  1  and  To  the  value  of  0.  These  values 
are  taken  as  a  matter  of  convenience.  It  is  not  the  specific  values 
that  are  of  importance  here.  The  emphasis  is  being  place  on  the 
dynamical  response  of  the  mathematical  functions. 

In  this  plot,  two  points  of  intersection  are  shown.  These  inter¬ 
section  points  correspond  to  the  possible  steady  state  conditions 


that  can  exist.  That  is,  the  states  where  the  rate  of  heat  generation 
equals  the  rate  of  heat  removal. 

In  passing  it  is  necessary  to  mention  that  three  possible  con¬ 
ditions  can  exist.  First,  the  slope  of  the  linear  function  can  be 
small  enough  such  that  the  linear  function  does  not  intersect 
the  exponential  curve.  In  that  case  no  steady  state  exists.  This 
means  that  the  system  response  would  rise  without  bound.  This 
situation  is  not  relevant  to  the  case  under  consideration.  Second, 
the  slope  of  the  linear  function  can  be  tangent  to  the  exponen¬ 
tial  curve.  This  is  actually  a  special  case  where  two  intersection 
points  converge  into  one  point.  The  third  case  is  the  situation 
shown  where  two  intersection  points  exist. 

3.3.  Discrete  dynamical  modeling 

To  analyze  this  situation,  an  approach  based  on  discrete 
dynamical  modeling  will  be  used  [7].  To  do  this,  the  steady 
state  expression  shall  be  rewritten  as 

In (*)  +  Tn  =  ln(m)  +  ln(r„+i  -To)  (1 1) 

or 

Tn+ 1  =  [k  exp(T^)  +  m7o]/m  (12) 

From  the  previously  assumed  values  of  m,  k  and  To,  the  following 
difference  equation  results: 

Tn+i  =  exp(Tn)/5  (13) 

It  is  now  possible  to  analyze  the  response  of  the  system  at  points 
other  than  the  two  intersection  or  equilibrium  points. 

The  analysis  begins  by  assigning  a  value  to  T  as  the  initial 
value.  This  value  determines  the  rate  of  heat  generation.  From 
this  assignment,  the  value  of  T  necessary  to  transfer  the  heat 
generated  to  the  environment  at  the  rate  at  which  it  was  produced 
is  determined.  Fig.  5  shows  the  consequence  of  assigning  an 
initial  value  of  T  that  is  between  the  two  intersection  points  (i.e., 
equilibrium  points).  This  diagram  is  sometimes  called  a  cobweb 
analysis.  The  consequence  of  this  calculation  is  that  the  T  is 
driven  downward  to  the  intersection  point  at  the  lower  value  of 
T.  The  figure  illustrates  this  effect. 


Fig.  5.  A  plot  of  the  linear  heat  transfer  process  and  the  exponential  Arrhenius 
function. 
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Fig.  6.  A  Cobweb  diagram  [7]  showing  the  sequential  steps  as  they  approach 
the  equilibrium  point. 

The  actual  iterative  calculation  is  shown  in  Fig.  6  showing 
the  convergence  to  the  equilibrium  point. 

The  convergence  take  place  at  about  T- 0.259.  The  same 
analysis  can  assume  that  the  initial  value  is  set  below  the  lower 
convergence  point.  If  one  assigns  the  initial  value  of  T  as  0.1, 
the  consequence  is  shown  in  Fig.  7.  This  time  the  rate  of  heat 
removal  obligates  the  temperature  for  heat  removal  to  rise  until 
the  convergence  point  is  reached  (Fig.  8). 

Finally,  an  initial  value  of  T  is  chosen  which  higher  than  the 
upper  equilibrium  point  the  result  is  shown  in  Fig.  9. 

The  starting  value  has  to  be  selected  near  and  above  the  equi¬ 
librium  point  to  see  that  T  diverges.  The  value  of  T  then  rises 
without  bound. 


Iteration  beginning  at  T=2.5 


Fig.  7.  Plot  of  the  calculated  points  as  the  system  approaches  the  equilibrium 
point  from  above. 


Iteration  beginning  at  T=0.1 


Fig.  8.  Plot  of  the  calculated  points  as  the  system  approaches  the  equilibrium 
point  from  below. 


Iteration  beginning  at  T=  2.555 


n 

Fig.  9.  Plot  of  the  calculated  points  as  the  system  deviates  starting  above  the 
upper  equilibrium  point. 

From  this  it  is  evident  that  the  two  equilibrium  points  have 
different  properties.  The  lower  valued  equilibrium  point  is  a 
fixed  point  attractor.  The  upper  point  is  often  called  a  saddle 
point  attractor.  If  one  can  achieve  the  exact  equilibrium  value 
of  T,  the  equilibrium  could  be  stable.  An  analogous  occurrence 
takes  place  at  the  two  equilibrium  points  of  a  pendulum.  The 
point  directly  below  the  fulcrum  is  a  fixed  point  attractor.  The 
equilibrium  point  directly  above  the  fulcrum  is  a  saddle  point 
attractor.  That  is,  it  is  an  unstable  equilibrium  point. 

Having  analyzed  the  behavior  of  the  functional  response  of 
the  controlling  equations,  an  interpretation  of  the  behavior  of 
thermal  runaway  can  be  given. 

3.4.  Discussion  of  thermal  runaway  using  this  modeling 
approach 

When  the  battery  is  being  charged  and  oxygen  evolution 
begins,  the  electrode  pack  begins  to  heat  up.  As  the  heating 
begins,  the  system  behaves  adiabatically.  Shortly  afterwards,  the 
heat  removal  process  begins.  The  rate  of  the  heat  removal  acts 
to  lower  the  temperature  within  the  electrode  pack.  From  the 
previous  discussion,  the  maximum  temperature  reached  within 
the  electrode  pack  becomes  very  important.  If  the  temperature 
never  exceeds  the  saddle  point  attractor  temperature,  the  system 
will  eventually  drop  to  the  lower  steady  state  equilibrium  point. 
However,  should  the  transient  temperature  rise  exceed  the  sad¬ 
dle  point  attractor  temperature,  the  heat  transfer  process  cannot 
remove  heat  faster  than  the  charging  current  can  generate  it.  A 
thermal  runaway  condition  exists. 

Without  attempting  to  quantify  the  actual  temperature  for  this 
cross  over  point,  it  is  instructive  to  consider  the  factors  that  can 
induce  a  high  initial  electrode  pack  temperature  and  the  factors 
affecting  the  curve  intersection  giving  the  saddle  point. 

3.5.  Factors  controlling  the  heat  capacity  of  the  electrode 
pack 

The  heat  capacities  of  the  constituent  parts  of  a  lead  acid 
battery  can  be  found  in  paper  by  Pavlov  [2].  The  weights  of 
the  electrode  pack  components  (the  active  block)  were  listed 
together  with  their  heat  capacities  of  those  components. 
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It  was  shown  that  the  electrolyte  solution  contributes  as  much 
as  85%  of  the  active  block  heat  capacity.  It  follows  that  small 
changes  in  the  electrolyte  quantity  can  significantly  affect  the 
heat  capacity  of  the  electrode  pack. 

Also,  it  necessarily  follows  that  the  loss  of  water  in  the  elec¬ 
trode  pack  can  affect  the  adiabatic  response  of  the  system.  By 
reducing  the  water  content,  an  applied  current  can  cause  an  exag¬ 
gerated  temperature  excursion  by  virtue  of  the  decrease  in  the 
heat  capacity  of  the  adiabatic  system. 

3.6.  Reducing  the  thermal  resistance  of  the  electrode  pack 

In  addition,  the  slope  of  the  curve  corresponding  to  the  rate  of 
heat  removal  depends  on  the  heat  transfer  coefficient  is  propor¬ 
tional  to  surface  area.  The  proportionality  constant  is  sometimes 
called  the  thermal  resistance.  Its  value  depends  on  the  insulative 
properties  of  the  material.  It  one  want  to  reduce  the  saddle  point 
intersection  point,  the  value  of  the  heat  transfer  coefficient  would 
be  reduced.  One  way  of  doing  that  is  to  increase  the  number  of 
dead  gas  spaces  within  the  electrode  pack.  That  is  most  easily 
accomplished  by  removing  water  from  the  electrode  pack. 

There  is  a  common  factor  affecting  the  thermal  runaway 
effect.  The  effect  of  water  loss  affects  thermal  runaway  by  virtue 
of  (a)  enhancing  the  transient  thermal  response  and  (b)  reducing 
the  steady  state  temperatures  where  the  rates  of  generation  and 
removal  are  equal. 

3. 7.  On  the  mechanism  of  the  thermal  runaway  process 

Based  on  the  preceding  discussion,  the  following  mechanistic 
interpretation  can  be  offered  (see  Fig.  10). 

During  a  normal  constant  voltage  charging  of  a  battery  as  one 
approaches  the  end  of  the  charging  process,  part  of  the  applied 
current  causes  oxygen  to  be  evolved.  Also,  during  overcharge, 
the  entire  current  is  converted  to  oxygen  production.  The  fol¬ 
lowing  mechanistic  steps  follow. 

The  water  is  electrolyzed  absorbing  storing  energy,  oxygen 
recombines  at  the  negative  electrodes  generating  heat  and  the 
Joule  heating  occurs  by  virtue  of  the  current  flow.  A  steady  state 
temperature  is  eventually  established  between  the  rates  of  heat 
generation  and  heat  removal.  Also,  a  steady  state  partial  pressure 


is  achieved  between  the  rates  of  oxygen  evolution  and  recombi¬ 
nation.  The  rise  in  temperature  results  in  an  autocatalysis  raising 
the  rate  oxygen  evolution  together  with  an  increased  cell  current 
flow  and  that  in  turn  produces  more  heat.  Also,  the  increased 
temperature  increases  the  rate  of  heat  removal. 

The  preceding  discussion  describes  the  situation  induced  by  a 
transient  thermal  response.  That  is,  should  the  initial  temperature 
rise  climb  to  level  where  the  rate  of  heat  removal  is  insufficient 
to  remove  the  heat,  then  the  autocatalysis  forces  the  reaction  to 
proceed  at  an  ever  increasing  rate. 

Of  course  this  process  cannot  raise  the  temperature  indefi¬ 
nitely.  Another  process  will  occur  to  limit  the  temperature  rise. 
One  possibility  is  that  the  excessive  oxygen  pressure  might  vent 
from  the  battery.  Another  possibility  is  that  the  electrolyte  might 
boil  and  remove  energy  via  a  liquid  to  gas  phase  transition. 

The  preceding  gives  a  picture  of  a  complex  mechanism  that 
allows  inferences  to  be  made.  These  inferences  permit  the  pro¬ 
posed  mechanistic  hypothesis  to  be  tested. 

At  this  point  in  the  exposition,  it  is  helpful  to  return  to  the 
observation  made  earlier.  Fig.  2  showed  an  interesting  attribute 
that  took  place  during  the  charging  process  when  the  battery  was 
in  thermal  runaway.  Over  a  period  of  time,  the  temperature  of  the 
battery  remained  essentially  constant  while  the  applied  current 
decreased.  This  constant  temperature  indicates  that  a  steady  state 
exists  between  the  rate  of  heat  generation  and  heat  removal.  The 
following  is  an  attempt  to  rationalize  this  in  view  of  the  observed 
decreasing  current.  To  do  this,  the  model  assuming  only  Joule 
heating  shall  be  invoked.  In  that  case  the  rate  of  heat  generation 
is  constant  so: 

r  R  =  k\  (steady  state  condition) 

In  this  case,  the  value  of  R  is  not  the  polarization  resistance.  It  is 
to  be  understood  as  the  effective  resistance  determined  as  if  the 
battery  were  dissipating  heat  as  an  ordinary  resistor.  A  resistor 
can  be  physically  analyzed  as 

R  =  pi/  A  =  fe/A  (at  constant  T)  (14) 

where  p  is  the  resistivity,  /  the  length  of  the  resistive  path  and 
A  is  the  cross-sectional  area.  The  resistivity  and  the  path  length 
are  held  constant  by  the  geometry  of  the  battery.  By  rearranging 


During  Constant  Voltage  Charging 
under  starved  electrolyte  conditions 
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Fig.  10.  A  representation  of  the  essential  steps  in  thermal  runaway. 
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and  lumping  the  constants  together 

I2k2/A  =  ki  (15) 

I2  =  k3A  (16) 

I  =  k47A  (17) 

The  inference  here  is  that  if  the  cross-sectional  area  of  the  resis¬ 
tor  changes  such  that  the  current  is  proportional  to  the  square 
root  of  the  cross-sectional  area,  a  constant  rate  of  heat  generation 
takes  place.  This  situation  can  conceivably  take  place  when  the 
electrolyte  in  the  inter-cell  gap  is  being  displaced  by  the  genera¬ 
tion  of  gas.  This  hypothesis  is  interesting  because  water  is  being 
moved  from  the  inter-cell  pack  but  not  from  the  cell  enclosure. 
This  provides  a  slightly  modified  interpretation  of  the  interpre¬ 
tation  given  by  Pavlov  [2].  Indeed,  the  low  water  loss  from  the 
battery  observed  as  a  consequence  of  the  thermal  runaway  can 
now  be  reconciled  with  the  changes  of  the  heat  capacity  of  the 
cell  pack  cause  by  water  loss.  That  is,  it  appears  that  electrolyte 
is  being  “pumped”  out  of  the  inter-electrode  gap  by  the  gas  evo¬ 
lution  process. 

The  mechanistic  hypothesis  is  adequate  for  further  testing. 
The  hypothesis  allows  for  designing  experiments  for  verifi¬ 
cation  and  provides  a  means  for  making  performance  predic¬ 
tions. 

3.8.  Carrying  the  analysis  further 

The  low  water  loss  and  the  presence  of  thermal  runaway 
suggest  that  the  increased  resistance  is  the  consequence  of  elec¬ 
trolyte  displacement  from  within  the  interelectrode  gap. 

Note  that  an  increased  interelectrode  resistance  caused  by 
electrolyte  movement  leads  to  an  interesting  speculation.  Elec¬ 
trolyte  movement: 

(1)  increases  the  IR  drop  in  the  inter-cell  gap; 

(2)  by  holding  constant  battery  voltage,  this  decreases  the  oxy¬ 
gen  electrode  potential; 

(3)  this  reduced  potential  decreases  the  rate  of  the  oxygen  gas 
evolution  reaction; 

(4)  this  in  turn  allows  the  electrolyte  to  flow  back  into  the  inter¬ 
electrode  gap; 

(5)  the  consequence  of  this  is  a  decrease  in  the  cell  resistance; 

(6)  the  rate  of  gas  evolution  increases  displacing  more  elec¬ 
trolyte; 

(7)  and  the  IR  drop  in  the  cell  increases  again. 

It  would  follow  that  by  maintaining  a  constant  applied  volt¬ 
age,  an  oscillating  current  is  a  possible  consequence.  This  effect 
has  been  previously  observed.  [Personal  communication  with 
John  Olson.] 

3.9.  Setting  up  dynamical  system  for  analysis 

From  the  preceding  description,  a  dynamical  model  can  be 
generated  by,  again,  studying  the  response  of  the  mathematical 


functions.  During  overcharge  the  applied  voltage  can  be  written 
as 

V(cell)  =  E( Pb)  +  E(  02)  +  IR  (18) 

where  E  is  the  electrode  potential.  Solving  for  the  oxygen  elec¬ 
trode  potential  gives 

E( 02)  =  y(cell)  -  E( Pb)  -  IR  (19) 

Since: 

(a)  V(cell)  is  a  constant  applied  voltage; 

(b)  Zs(Pb)  is  depolarized  by  O2  recombination; 

then 

E( 02)  =  kx  -  IR  (20) 

and 

rate  of  reaction  =  k 2  exp  k^{E  (O2))  (21) 

Substituting 

rate  of  reaction  =  k 2  exp  k/±(—IR)  (22) 

Next,  let  IR  =  pT ,  where  p  shall  be  considered  a  coupling  con¬ 
stant.  The  general  properties  of  this  coupling  are  that  p  is  small 
in  an  environment  where  the  electrolyte  easily  flows  back  into 
replace  the  gas  bubbles  produced  such  as  in  a  flooded  system  and 
large  in  the  case  where  it  is  difficult  to  displace  the  gas  bubbles 
such  as  in  a  immobilized  electrolyte  condition.  Then 

I  (rate  of  reaction)  =  k 2  exp  k/±(—pT)  (oxygen  electrode) 

(23) 
and 

I  (rate  of  reaction)  =  £5  exp  ke  T  (Arrhenius  component) 

(24) 

These  equations  illustrate  the  competing  processes  of  the 
accelerating  effect  of  the  Arrhenius  response  and  the  retard¬ 
ing  effect  of  the  decreased  polarization  caused  by  the  IR 
drop. 

By  setting  up  the  analysis  graphically  as  a  discrete  dynamical 
model,  some  inferences  can  be  made  about  the  system  behavior 
based  on  this  model.  The  first  case  shall  assume  that  a  weak  cou¬ 
pling  takes  place  between  the  thermally  catalyzed  gas  production 
and  the  ability  of  the  electrolyte  to  return  into  the  inter-cell  gap 
to  reduce  the  IR  drop.  The  IR  drop  reduction  serves  to  reduce 
the  retarding  influence  on  the  electrode  potential  driven  oxygen 
evolution  reaction. 

The  abscissa  (the  v-axis)  is  proportional  to  the  temperature 
and  the  ordinate  (y-axis)  is  proportional  to  the  rate  of  the  reaction 
(or  applied  current).  The  iteration  represents  initially  the  effect 
of  temperature  on  the  temperature  catalyzed  reaction  rate.  To 
sustain  that  rate,  the  IR  drop  required  to  do  this  would  have  to 
rise.  That  is,  IR  is  determined  by  the  coupling  constant.  At  this 
temperature,  the  thermally  catalyzed  process  would  increase  the 
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Fig.  11.  x-Axis  is  a  function  of  T  and  the  _y-axis  is  a  measure  of  the  reaction  rate  p  ~  weak  coupling. 


rate  of  reaction  (i.e.,  the  gas  generation).  To  sustain  that  rate,  the 
IR  drop  would  have  to  drop  corresponding  to  a  lower  temper¬ 
ature.  This  repetitive  sequence  of  events  is  seen  to  converge  to 
an  equilibrium  point  (an  attractor).  Fig.  11  illustrates  that  the 
equilibrium  point  is  stable  and  the  convergence  can  take  place 
even  more  quickly  as  the  coupling  becomes  weaker. 

In  the  case  of  a  strong  coupling  condition,  the  equilibrium 
point  is  unstable  (a  saddle  point  attractor).  Fig.  12  shows  that 
the  temperature  excursions  exhibit  a  reinforced  oscillation  that 
grows  without  bound.  Clearly,  this  suggests  that  some  other  pro¬ 
cess  would  intervene  so  as  to  eventually  limit  the  temperature 
growth. 

Fig.  13  shows  an  intermediate  condition  which  is  the  transi¬ 
tion  between  a  stable  and  unstable  equilibrium  point. 

A  “limit  cycle”  appears  around  the  equilibrium  point.  This 
represents  a  special  condition  where  the  system  can  develop  a 
sustained  oscillation  that  neither  move  toward  or  away  from  the 
equilibrium  point. 


3.10.  Experimental  observations 

These  model  predictions  provided  the  guidance  for  search¬ 
ing  for  the  effects  experimentally.  The  chosen  experimental 
allows  the  prediction  of  this  model  to  be  tested.  The  methodol¬ 
ogy  for  inducing  the  thermal  runaway  experimentally  permitted 
the  number  of  constant  current  cycles  used  to  induce  the  effect 
to  be  altered.  Fifty  cycles  was  found  to  be  sufficient  to  induce 
the  most  severe  effect.  To  establish  a  mild  pretreatment,  eight 
constant  current  cycles  was  applied  to  a  new  battery  that  was 
followed  by  the  usual  constant  voltage  charge.  This  resulted  in 
the  data  shown  in  Figs.  14  and  15.  This  represents  a  condition 
corresponding  to  a  loose  coupling.  The  plots  that  follow  include 
the  response  as  a  function  of  time  and  accumulated  charge.  This 
is  the  same  data  represented  in  two  different  ways. 

The  system  appears  to  have  stabilized  showing  no  apparent 
transition  to  a  thermal  runaway.  A  careful  look  at  the  data  sug¬ 
gests  a  possible  onset  but  the  stabilization  is  evident. 
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Fig.  12.  x-Axis  is  a  function  of  T  and  the  _y-axis  is  a  measure  of  the  reaction  rate  p  ~  strong  coupling. 
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Fig.  13.  x-Axis  is  a  function  of  T  and  the  y-axis  is  a  measure  of  the  reaction  rate  p  ~  medium  coupling. 


CV  Charging 

Mild  Pretreatment 


Fig.  14.  Consequence  of  loose  coupling  shown  as  a  function  of  time. 


An  additional  eight  constant  current  cycles  was  applied  to 
this  same  battery  followed  by  the  usual  constant  voltage  charge. 
In  this  case,  the  effect  of  oscillation  was  observed.  The  data 
appears  in  Figs.  16  and  17. 

The  observed  oscillatory  behavior  appears  to  be  chaotic.  It  is 
consistent  with  the  notion  that  it  occurs  when  a  special  situation 
takes  place  consistent  with  an  intermediate  coupling  between 
the  electrolyte  displacement  by  gas  evolution  and  the  return  to 
the  electrolyte  to  reduce  the  IR  drop  in  the  inter-cell  gap. 


CV  Charging 

Mild  Pretreatment 


Fig.  15.  Consequence  of  loose  coupling  shown  as  a  function  of  charge. 


Constant  Voltage  Charging 


Fig.  16.  Consequence  of  loose  coupling  shown  as  a  function  of  time. 


Finally,  a  full  50  cycles  was  applied  to  drive  the  system  into  its 
worse  case  condition.  This  would  represent  the  situation  where  a 
tight  coupling  would  exist.  The  data  is  shown  in  Figs.  18  and  19. 

It  is  to  be  noted  that  the  predicted  reinforced  oscillations  did 
not  occur.  Since  the  uncontrolled  oscillation  represents  the  pre¬ 
dicted  behavior  and  it  did  not  occur,  the  battery  was  examined 
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Fig.  17.  Consequence  of  loose  coupling  shown  as  a  function  of  accumulated 
charge. 
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Fig.  18.  Consequence  of  tight  coupling  shown  as  a  function  of  time. 


CV  Charging  During  TR 
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Fig.  19.  Consequence  of  tight  coupling  shown  as  a  function  of  accumulated 
charge. 

closely  to  attempt  an  assessment  as  to  why  the  ever  increas¬ 
ing  oscillations  were  absent.  A  close  examination  showed  that 
all  six  vents  on  the  battery  were  open  with  steam  passing  from 
the  battery  into  the  environment.  As  was  mentioned  earlier,  the 
temperature  measurement  shown  was  the  temperature  of  the  neg¬ 
ative  electrode  terminal  with  the  heavy  connectors  and  cabling 
attached.  The  temperature  inside  the  battery  is  necessarily  higher 
than  the  terminal  temperature.  The  presence  of  the  steam  sug¬ 
gests  that  localized  electrolyte  boiling  was  taking  place.  Also, 
a  very  strong  stench  of  hydrogen  sulfide  was  being  emitted 
from  the  battery.  From  these  observations  one  can  conclude  that 
secondary  processes  take  place  acting  to  prevent  a  further  rise 
in  the  battery  temperature.  The  predicted  reinforced  oscillation 
appeared  to  have  been  damped  by  secondary  processes. 

The  generation  of  hydrogen  sulfide  requires  some  additional 
remarks.  The  thermodynamic  situation  in  a  lead-acid  battery 
indicates  that  the  reduction  of  sulfuric  acid  by  lead  is  sponta¬ 
neous.  However,  the  kinetics  of  that  process  are  very  slow  under 


normal  operating  conditions.  Further,  any  trace  amounts  pro¬ 
duced  can  be  effectively  scrubbed  by  the  oxidizing  power  of 
the  positive  electrode.  In  the  experiments  performed,  the  high 
temperature  produced  in  the  interelectrode  gap  was  sufficient  to 
catalyze  this  kinetically  slow  process  to  a  level  where  the  hydro¬ 
gen  sulfide  generation  was  very  apparent  the  people  standing  in 
the  immediate  vicinity. 

4.  Conclusions 

The  cause  of  the  thermal  runaway  appears  to  be  associated 
with  processes  controlled  by  the  separator.  The  approach  to  mini¬ 
mizing  the  TR  effect  would  be  to  allow  for  a  more  efficient  return 
of  the  displaced  electrolyte  into  the  inter  cell  gap.  There  is  also  a 
potential  here  for  using  proton  conductive  coatings  on  the  fiber 
glass  surface. 

Temperature  control  could  be  used  to  place  limits  on  the  cur¬ 
rent  allowed  to  pass  into  the  cell.  This  would  allow  more  time 
to  permit  electrolyte  to  reenter  the  inter-cell  gap. 

From  the  model  presented  above,  it  would  follow  that  hydro¬ 
gen  generation  from  a  negative  limited  charging  condition  could 
generate  a  thermal  runaway.  That  is,  the  electrolyte  displacement 
can  be  brought  about  by  any  gas  generation  process. 

It  is  also  noteworthy  that  the  mechanism  is  independent  of 
cell  chemistry  and  should  apply  to  cells  having  a  constrained 
electrolyte  that  can  be  displaced  by  gas  generated  during  charg¬ 
ing. 

Also,  there  exists  a  problem  associated  with  the  catalysis 
of  other  passive  processes.  H2S  production  and  in  Li-ion  cells, 
electrolyte  conflagration  are  examples.  The  TR  effect  can  serve 
to  trigger  an  otherwise  kinetically  slow  reaction. 
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